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Abstract. We report detailed follow-up observations of the cataclysmic variable HS 2331+3905, identified as an 
emission-line object in the Hamburg Quasar Survey. An orbital period of 81.08 min is unambiguously determined 
from the detection of eclipses in the light curves of HS 2331+3905. A second photometric period is consistently 
detected at P ~ 83.38 min, ~ 2.8% longer than P or b, which we tentatively relate to the presence of permanent 
superhumps. High time resolution photometry exhibits short-timescale variability on time scales of ~ 5 — 6 min 
which we interpret as non-radial white dwarf pulsations, as well as a coherent signal at 1.12 min, which is likely to 
be the white dwarf spin period. A large-amplitude quasi-sinusoidal radial velocity modulation of the Balmer and 
Helium lines with a period ~ 3.5 h is persistently detected throughout three seasons of time-resolved spectroscopy. 
However, this spectroscopic period, which is in no way related to the orbital period, is not strictly coherent but 
drifts in period and/or phase on time scales of a few days. Modeling the far-ultraviolet to infrared spectral 
energy distribution of HS 2331+3905, we determine a white dwarf temperature of Teff — 10 500 K (assuming 
A/ W d = 0.6 Mq), close to the ZZ Ceti instability strip of single white dwarfs. The spectral model implies a distance 
of d = 90 ± 15 pc, and a low value for the distance is supported by the large proper motion of the system, 
/i = 0.14" yr~ . The non-detection of molecular bands and the low J, H, and K fluxes of HS 2331+3905 make 
this object a very likely candidate for a brown-dwarf donor. 
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1. Introduction 

In the current stand ard model of cataclysmic variable 
(CV) evolution (see iKind fl988^ the systems evolve to 
shorter orbital periods, reach a minimum period at ~ 
70 min and then begin to evolve back to longer orbital 
periods. As a CV approaches the minimum period, the 
mass of the donor star decreases until it can no longer sus- 
tain hydrogen burning. The low mass secondary star be- 
comes i ncreasingly degenerate, entering the brown dwarf 
regime. iKblbl l)l993|) estimated that 70% of all CVs should 
have such "substellar" secondaries. Despite these large 
predicted numbers, only a handful of CVs are currently 
believed to harbour a brown dwarf donor, with WZ Sge, 
EG Cnc and EF Eri being among the best candidates 
llLittlefair et alJl2003t IPattersonl Il99i iBeuermann et alJ 
2000). So far, for none of them does a compelling proof of 
the degenerate nature of the star exist. The discrepancy 
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between the observed and predicted numbers of brown 
dwarf CVs could be due to selection effects: because of the 
low mass transfer rates predicted for these systems they 
may be substantially fainter than most CVs and may have 
very rare outbursts. Alternatively, it may also be that the 
standard theory is wrong , and that CVs have either not 
reached the period minimum wher e their donors become 
degenerate ijKing fc Schenkerll2002f) . o r they may even die 
before reaching the period minimum l|Pattersonll 1 99Sj) . 

In order to be able to test properly any of these 
hypotheses we need to reduce the observational bi- 
ases that afflict the current CV population to a min- 
imum. This is the main purpose of an undergoing 
large-scale sea rch for CVs using the Hamburg Quasar 
Survey fflQS: iHageu et al.l Il995t ICausicke et all 2000, 
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200 21 ICansicke et al.ll2004t lArauio-Betancor et all I2OO.I ; 
Rodriguez-Cil et al.ll2004|) . One of our most recent discov- 
eries, HS 2331+3905 (HS 2331 thereafter), was selected as 
a CV candidate on the basis of its spectral characteristics 
in the HQS. Specifically, its HQS spectrum displays the 
strong Balmer emission lines that suggest ongoing accre- 
tion in CVs. 

In this paper we present the results of follow-up spec- 
troscopy and photometry of HS2331 obtained over three 
years, which confirm the CV nature of this object and 
unveil a very unusual picture. In Sec.[21 we describe the 
photometric and spectroscopic (ground and space) obser- 
vations and data reduction and calibration. The analysis 
of the photometry and ground-based spectroscopy is pre- 
sented in Sec. 13 and Q] In Sec.0] we also present a model 
for the spectral energy distribution of HS 2331, from which 
we derive a temperature for the white dwarf and distance 
to the system. Finally, in Sec. 03 we use our findings in 
trying to construct a preliminary picture of the nature of 
HS2331. 

2. Observations and Data Reduction 

2.1. Photometry 

2.1.1. Braeside Observatory 

Differential i?-band CCD photometry of HS 2331 was ob- 
tained at the Braeside Observatory immediately following 
the identification of its CV nature in September-October 
2000 fSect. l2~2~T|l . using a 0.4 m reflector equipped with 
a SITe512 CCD camera (Table[TJ). HS2331 was again ob- 
served from Braeside in November 2001 and October 2003, 
on both occasions in white light. The data were bias- 
subtracted, dark current-subtracted and flat-fielded in a 
standard fashion, and instrumental aperture magnitudes 
of HS 2331 were derived relative to USNO-A2.0 1275- 
18486676 {R = 14.6, labelled 'CI' in Fig.[TJ. The mean 
magnitude of HS 2331 during the 2000 observations was 
R ~ 16.5. 

The Braeside CCD light curves of HS 2331 display peri- 
odic variability with an amplitude of ~ ±0.1 mag (Fig.^J, 
which has consistently been detected in all three years cov- 



Fig. 1. 10' x 10' Finding chart of HS 2331, obtained from 
the Digitized Sky Survey 2. The coordinates of the CV 
are a(2000) = 23^34 m 01.6 s , (5(2000) = +39°21'41.4". 
Differential magnitudes of HS 2331 were obtained relative 
to the comparison stars CI (USNO-A2.0 1275-18486676) 
and C2 (CSC 0323100595). 



ered by our observations. The morphology of the photo- 
metric modulation is best described by a double-humped 
pattern with a period of ~ 80 min. The relative strength of 
the two humps varies substantially between the individual 
nights. The 2003 data reveal narrow absorption dips cen- 
tred on some of the observed minima between the humps. 

2.1.2. Kryoneri Observatory 

Filterless CCD photometry of HS 2331 was obtained dur- 
ing two nights in October 2002 and during 7 nights in 
August 2003 at the 1.2 m Kryoneri telescope using a 
SI-502 516 x 516 CCD camera (TableHJ. The Kryoneri 
data were reduced fo llowing the procedure described in 
iGansicke et alJ ((2004) , and differential magnitudes were 
derived relative to CI (Fig.GJ. 

The October 2002 light curves, obtained with a time 
resolution of 25 s - 30 s, were essentially similar to 
the Braeside observations (Fig.|5J), displaying a double- 
humped structure with a period of ~ 80 min. The higher 
time resolution data obtained in August 2003 provide 
stronger evidence for the absorption dips detected in the 
Braeside observations, and clearly reveal additional short- 
term variability on a time scale of ~ 5 min (e.g. August 
14; Fig.|5J. In addition the changes of the shape of the 
double-humped structure are more pronounced than in the 
Braeside observations. On August 19, one of the humps is 
extremely weak, on August 20 both humps have nearly 
equal amplitude. 



2.1.3. ULTRACAM on the William Herschel Telescope 

The highest time resolution and signal-to-noise ratio 
CCD photometry was obtained on November 10, 2003 
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Photometry 

Date UT Time Filter Exp. [s] # Frames 



40 cm Braeside Observatory 



2000 Sep 25 


04:35 - 


12:29 


R 


50 


510 


2000 Sep 28 


05:26 - 


10:54 


R 


85 


211 


2000 Oct 03 


02:21 - 


12:01 


R 


95 


338 


2001 Nov 21 


04:13 - 


08:44 


C 


85 


152 


2001 Nov 22 


02:43 - 


08:45 


C 


55 


302 


2003 Oct 12 


02:19 - 


11:44 


C 


55 


552 


2003 Oct 13 


02:36 - 


11:30 


c 


55 


542 


2003 Oct 15 


03:56 - 


11:23 


c 


55 


446 



1.2 m Kryoneri Observatory 



2002 Oct 15 


18:41 - 


22:39 


c 


30 


340 


2002 Oct 16 


17:29 - 


22:13 


c 


25 


429 


2003 Aug 14 


22:56 - 


02:52 


c 


10 


1029 


2003 Aug 15 


23:14 - 


02:57 


c 


10 


1012 


2003 Aug 16 


23:33 - 


03:06 


c 


10 


980 


2003 Aug 17 


23:32 - 


03:04 


c 


10 


948 


2003 Aug 18 


22:23 - 


03:01 


c 


10 


1249 


2003 Aug 19 


22:09 - 


03:04 


c 


10 


1311 


2003 Aug 20 


21:56 - 


03:07 


c 


10 


1432 



4 m WHT/ULTRACAM 

2003 Nov 10 19:23 - 22:00 u 1 10013 

2003 Nov 10 19:23 - 22:00 g' 1 10013 

2003 Nov 10 19:23 - 22:00 r' 1 10013 



Spectroscopy 

Date UT Time Grating Exp. [s] # Spectra 



2.2 m Calar Alto/CAFOS 



2000 Sep 20 


00:50 


B-200 


600 


1 


2000 Sep 20 


01:05 


R-200 


600 


1 


2000 Sep 23 


21:20 - 02:30 


B-100 


Ron 


9 1 


2000 Sep 24 


00:22 - 03:37 


B-100 


600 


16 


2000 Sep 24 


03:47 - 04:27 


R-100 


600 


3 


2003 Sep 11 


22:16 - 04:20 


G-100 


600 


28 


2003 Sep 15 


02:28 - 04:35 


G-100 


600 


12 


2003 Sep 16 


23:48 - 04:04 


G-100 


600 


23 


2.5 m Isaac Newton Telescope/IDS 






2002 Aug 26 


02:37 - 03:60 


R632V 


600 


9 


2002 Aug 28 


00:33 - 01:56 


R632V 


600 


9 


2002 Sep 01 


04:16 - 04:47 


R632V 


600 


4 


2002 Sep 02 


03:09 - 03:41 


R632V 


600 


4 


Hubble Space Telescope/STIS 






2002 Oct 24 


11:34 


G140L 


700 


1 


2.4 m MDM Hiltner 








2003 Oct 12 


07:28 - 08:15 


600 


420 


7 


2003 Oct 13 


07:38 - 08:17 


600 


420 


6 


2003 Oct 14 


08:03 - 08:59 


600 


420 


8 


2003 Oct 15 


02:51 - 10:48 


600 


420 


21 


2003 Oct 16 


03:45 - 04:40 


600 


420 
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using the 3-beam ULTRACAM CCD camera (which 
takes CCD data in three colour channels simultaneously, 
iDhillon et afll2002|) on the William Herschel Telescope. 
HS2331 was observed with a Is time resolution, using 
Sloan u', g' , and r' filters. The data were reduced with the 
ULTRACAM pipeline. Differential magnitudes of HS2331 
were obtained relative to GSC 0323100595 (labelled 'C2' 
in Fig.^l, and converted to apparent magnitudes using 
observations of the Sloan standard SA100-280, which was 
observed at the end of the night. The mean magnitudes 
of HS2331 were u' = 16.1, g' = 16.3 and r' = 16.0. These 
values are subject to some small systematic uncertainty 
related to the yet uncalibrated colour terms accounting 
for the differences between the ULTRACAM and Sloan 
setups. 

The ULTRACAM light curves impressively confirm 
the presence of variability on a time scale of ~ 5 min seen 
in the Kryoneri data, as well as additional variability on 
even shorter time scales (Fig.^J. 

2.2. Optical Spectroscopy 
2.2.1. Calar Alto 

HS 2331 was identified as a CV from low-resolution spec- 
troscopy obtained at the Calar Alto 2.2 m telescope on 
September 20 2000. The CAFOS focal reductor spectro- 
graph was used in conjunction with the standard SITe 
CCD. A single pair of blue/red identification spectra was 



obtained with the B-200 and R-200 grating through a 2" 
slit (see TableQ]). The wavelength range of the combined 
spectra extends from 3500 to 10000 A at a spectral reso- 
lution of approximately 10 A. The online reduction of the 
identification spectra using the CAFOS quick look con- 
text within MIDAS immediately revealed the CV nature 
of HS 2331 by the presence of broad double peaked Balmer 
and He I emission lines in the spectra. 

Time- resolved follow-up spectroscopy of HS 2331 was 
carried out during the same observing run, on September 
23 and 24 2000. A total of 43 spectra of 600 s each 
were obtained during the two nights, which were mainly 
photometric except for very thin cirrus during the last 
~ 1.5 hours of the second night (TableQ]). The B-100 grism 
was used along with a 1.5" slit to obtain 40 spectra with 
a useful wavelength range of 3500-6300 A and a spectral 
resolution of ~ 4 A. In addition, 3 red spectra were ac- 
quired during the second night, using the R-100 grism 
with a 1.5" slit width. The red spectra cover a wave- 
length range of 6000 to 9200 A at a spectral resolution 
of ~ 4 A. Mercury-cadmium, helium-argon and rubidium 
arcs were taken regularly throughout both nights to pro- 
vide for the wavelength calibration. Data from the flux 
standards G191B2B and £Cas were also obtained to cor- 
rect for the instrument response of the "blue" and "red" 
spectra respectively. The CAFOS acquisition image were 
taken with a Johnson V filter in order to correct the spec- 
tra for slit losses. 




HJD 



2927.9 
2 450 000 



2872.6 
HJD 



2872.55 
2 450 000 



Fig. 2. Sample light curves of HS2331 based on differential CCD photometry obtained at the 40 cm Braeside 
Observatory (left) and at the 1.2m Kryoneri Observatory (right). 
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Fig. 3. Light curve of HS 2331 obtained at a time resolution of 1 s with the ^'-filter using ULTRACAM on the WHT. 



All spectra were reduced in a standard manner us- 
ing the Figaro package within the Starlink software col- 
lection. The bias level was subtracted from the calibra- 
tion and object images using the median from a set of 
unexposed frames. Dome fiat field images, taken with 
the same set up as the object frames, were used to re- 
move pixel to pixe l variations. The spectra were then op- 
timally extracted l|Hornelll986j) and sky line subtracted. 
A fourth-order polynomial was fitted to obtain the dis- 
persion relation for each of the arc spectra, giving an 
RMS smaller than 0.1 A in all cases. Each target spec- 
trum was then wavelength-calibrated by interpolating two 
neighboring arc spectra to account for the possible drift 
of the wavelength-to-pixel scale. The instrumental re- 
sponse function was computed as the ratio of a spline 
fit to the line-free continuum of the B-100/R-100 spec- 
tra of the flux standards to the flux- calib r ated d ata from 
iFilippenko fc Greensteinl l)l984|) andlOkd jl990). All tar- 



get spectra were flux-calibrated and corrected for atmo- 
spheric extinction using this response function. 

We used the acquisition images obtained in con- 
junction with the CAFOS spectroscopy to estimate the 
brightness of HS 2331 during our Calar Alto observa- 
tions. Aperture photometry relative to CSC 0323100595 
(V = 12.75) located ~ 35" south-east of HS 2331 ('C2' in 
Fig.QJ resulted in V ~ 16.4, which is consistent with the 
magnitude obtained from the HQS direct plate, and with 
the brightness level of all our photometry fSect. l2TT| l. A 
grand average spectrum of HS 2331 was created by com- 
bining the mean of all "blue" and "red" spectra from the 
flux-calibrated (and corrected for atmospheric extinction) 
B-100/R-100 data set. This average spectrum was then 
scaled to V — 16.4, as determined from the CAFOS ac- 
quisition images to account for possible slit losses (Fig.^J . 
The applied correction of ~ 0.2 mag is comparable to the 
intrinsic variability of HS 2331, and indicates that the ab- 
solute error in our flux calibration is < 20%. 
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Table 2. EW and FWHM (corrected for the instrumental 
resolution) of the strongest Balmer emission lines and of 
the strongest FUV lines. 



Line 


EW [A] 


FWHM [km/s] 


FWHM [A] 


Ha 


185 ± 14 


1312 ± 5 


28.7 ±0.1 


H/3 


75 ±9 


1475 ± 6 


23.9 ±0.1 


H>, 


40 ±4 


1583 ± 7 


22.9 ± 0.1 


H<5 


30 ±13 


1403 ± 7 


19.2 ± 1.1 


Civ 


234 ± 14 


2478 ± 19 


12.8 ± 0.1 


Si iv 


52 ± 13 


3749 ± 236 


17.5 ±1.1 


Cn 


67 ± 12 


2476 ± 68 


11.0 ±0.3 


Cm 


111 ± 19 


2500 ± 51 


9.8 ±0.2 


Nv 


42 ± 15 


2417 ± 48 


10.0 ±0.2 



A further set of observations were obtained in 
September 2003 at the 2.2 m Calar Alto telescope. This 
time the G-100 grism was used which together with a slit 
width of 1.2" provided a wavelength range of 4240-8300 A 
and a spectral resolution of ~ 4.1 A (Tablc^J. These spec- 
tra were reduced in an analogous fashion as described 
above. 

The spectrum of HS2331 (Fig.gJ displays broad 
double-peaked Balmer and Hei emission lines, suggesting 
an origin in the accretion disc. The higher series of the 
Balmer lines are flanked by extremely broad absorption 
troughs, which are reminiscent of the Stark-broadened ab- 
sorption lines observed in the high-gravity atmospheres of 
cool white dwarfs (the width of the absorption lines are 
too wide to be of accretion disc origin). The red part of 
the spectrum does not contain any spectral features that 
could be ascribed to the emission of the secondary star. 

TableEJlists the equivalent widths (EW) and full width 
at half maximum (FWHM, corrected for the instrumen- 
tal response), of Ha-EW. The EW and the FWHM mea- 
surements were obtained by fitting a single Gaussian in 
all cases. Average spectra were used for the Balmer lines: 
Calar Alto 2000 for H<5 and Calar 2003 for Ha-H<5. 



2.2.2. 2.5 m Isaac Newton Telescope 

A total of 26 spectra of 600 s each, spread out over a week, 
were obtained in August-September 2002 (TableHJ. We 
used the Intermediate Dispersion Spectrograph (IDS) on 
the 2.5m Isaac Newton Telescope (INT), at the Roche 
de los Muchachos observatory in La Palma. The IDS was 
equipped with the R632V grating and the 2048 x 4100 
pixel EEV10 detector. Using a slit width of 1.5" the setup 
provided an unvignetted wavelength range of - 4400 A - 
6800 A and a spectral resolution of ~ 2.3 A. Copper- Argon 
wavelength calibrations (arcs) were obtained at the begin- 
ning and end of each observation block of HS2331. The 
reduction procedure applied to this data set is identical to 
the one applied to the Calar Alto/CAFOS data above. 



2.2.3. 2.4 m Hiltner telescope 

Most recently, we observed HS 2331 from the 2.4 m Hiltner 
telescope at MDM observatory (Tabled . We obtained a 
total of 50 spectra of 420 s each in October 2003. The 
modular spectrograph combined with the 600 line mm -1 
grating and a SITe 2048 x 2048 CCD detector, yielded a 
spectral dispersion of ~2.0 A pixel -1 from 4000 to 7500 A. 
The reducti on technique applied to t his set of data is as 
explained in IThorstensen et al.l 1 1998j) . 



2.2.4. Hubble Space Telescope 

Far-ultraviolet (FUV) spectroscopy of HS 2331 was ob- 
tained with the Hubble Space Telescope/Space Telescope 
Imaging Spectrograph (HST/STIS) on October 24, 2002 
as part of an ongoing snapshot survey of CVs. The ob- 
servations were obtained using the G140L grating and 
the 52" x 0.2" aperture, providing a spectral resolution 
of R w 1000 over the wavelength range 1150 - 1710A. 
The data were pipeline-processed within IRAF using 
CALSTIS V2.13b. The F28x50LP magnitude of HS2331 
measured from the STIS CCD acquistion image confirms 
that the system was at its usual brightness level dur- 
ing the FUV observations. The STIS spectrum (Fig.EJ 
contains the mixture of low and high-ionisation emis- 
sion lines often observed in CVs: Cm A 1176, Sim A 1206, 
NvAA 1239,43, CnA1335, a broad complex of Sim near 
1300 A, Si iv AA 1394,1403, CivAA 1548,51 and probably 
He ii A 1640. The emission line flu x ratios are within 
the normal range ob served in CVs ijMauche et al.l [l"997t 
iGansicke et al.ll2003j) . Several of the lines contained in the 
STIS are double-peaked, again suggesting an origin in an 
accretion disc. The EW and FWHM of the strongest lines 
in the STIS spectrum were estimated by fitting a Gaussian 
to the observed line profiles, and the results are reported 
in Tabled Correcting the FWHM for the separation of the 
Cn, Civ and Siiv doublets results in line widths which 
are comparable to those of the Balmer lines, suggesting a 
common origin of the FUV and optical emission lines. 

The continuum flux underlying the emission lines is 
slightly increasing towards the red, with a sudden upturn 
in flux at wavelengths A <; 1650 A. 

3. Analysis: Photometry 

We have analysed the complex variability detected in the 
light curves of HS2331 (for a small sam ple see F ig. El and 
13) by computing Scargle periodograms (Scarglc 1982) for 
the entire CCD photometry of HS2331, as well as for in- 
dividual sub-sets of the data, using the MIDAS/TSA con- 
text. Prior to this period analysis, we have subtracted the 
nightly mean of each light curve to account for night-to- 
night variations in the overall brightness of the system and 
more importantly of the different detectors/filters used 
during the observations. The periodogram of HS 2331 cal- 
culated from the Kryonery 2003 data, the largest and best 
sampled photometric data set, (Fig.® top panel) contains 
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Fig. 4. Mean, flux-calibrated "blue" and "red" spectra Fig. 5. HST/STIS spectrum of HS 2331. 
(Calar Alto 2000, B-100/R-100 data set) of HS 2331 com- 
bined together. The arrows indicate the wavelength range 
covered by the B-100 and R-100 gratings. 



a multitude of signals, as anticipated by our visual inspec- 
tion of the light curves fSect. l2.ljl . 

The strongest peak is located at 35.52 d _1 (40.54 min), 
a weaker signal is centred on ^17.27 d _1 (83.38 min), and 
some power is also found at 53.28 d _1 (27.03 min) and 
71.04d _1 (20. 27 min). In the high-frequency range there 
is a strong peak at 1284. 7d _1 (1.12min) and a weaker one 
at 256.6 d^ 1 (5.61 min). 

3.1. The orbital period 

Fig.EI shows a close-up view of the Scargle periodogram 
of the ~40min signal (middle left panel). The strongest 
peak is located at 40.542632 ± 0.000008 min (35.518168 ± 
0.000007 d _1 ) with the one a error obtained from a sine 
fit to the entire photometric data. The other peaks seen 
in the periodogram are aliases separated by multiples of 
Id -1 . The sharpness and power of the 40.54 min signal 
indicates a very high degree of coherence, suggesting that 
it is related to a stable "clock" within the binary system, 
such as the binary orbital period or the white dwarf spin 
period. 

Folding the entire 4 years worth of photometry over 
exactly twice the 40.54min signal, i.e. 81.08 min, results 
in an average light curve (Fig.UJ) that unambiguously 
recreates the morphology observed in the individual light 
curves (Fig.[21 & a double- humped modulation with 
a narrow dip. Dips of such shape are observed either in 
strongly magnetic CVs, when the magnetically funnelled 
accretion stream passes in front of the white dwarf, or as 
grazing eclipses of the hot spot in the accretion disc in 
non-magnetic CVs. Stream eclipses in magnetic CVs are 
relatively unstable features as the location of the stream 
within the binary frame changes as a function of the mass 
transfer rate and therefore we favour the latter scenario 
as the cause of the narrow dips. 



The light curve of HS2331 (Fig.EI&CJ bears a striking 
resemblance to that of one of the most (in)famous and 
thoroug hly studied short per iod CVs: WZSge (see e.g. 
Fig. 1 of IPatterson et al . 1998). Both objects exhibit dou- 
ble humps in their light curves, in both cases one of the 
humps appear to suffer from sporadic absorption events, 
and both systems show a sharp and relatively shallow 
eclipse. From the strict coherence of the 40.54 min, the 
detection of periodically recurrent eclipses, and the close 
analogy to WZ Sge, we conclude that the orbital period of 
HS2331 is 81.08 min, and we define the following eclipse 
ephemeris 

0o = HJD 2451812.67765(35) + 0.05630921(1) x E (1) 

where phase zero corresponds to the time of mid eclipse. 
Higher harmonics of the orbital period are clearly visible 
in the periodogram (Fig.[fjJ) at 27.03 min (53.28 d -1 ) and 
20.27min (71.04d~ 1 ). 

3.2. Quiescent superhumps? 

Only a weak signal is detected at the orbital period of 
81.08 min (Fig[S] middle right). This is not too much of 
a surprise because of the double-humped structure of the 
orbital modulation, which shifts most power to half the or- 
bital period. Surprising is, however, the detection of a rel- 
atively strong signal close to P ar b, with a period of either 
~ 83.38 min (~ 17.27 d" 1 ) or ~ 88.51min (~ 16.27d~ 1 ). 
Both these 1-day aliases are of similar strength, and their 
relative power varies between different observing runs. In 
contrast to the 35.52 d _1 orbital signal (Fig.Elmiddle left), 
the periodogram around the 83.38 min/88. 51 min signals is 
extremely complex, displaying significant substructure in 
excess to the window function, suggesting either the pres- 
ence of more than one periodicity and/or a lower degree 
of coherence. Figure|Hl shows the light curve folded over 
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the 83.38 min signal detected in the periodogram of Fig.© 
(middle- left panel). The absence of eclipses in this folded 
light curve clearly rules out the 83.38 min signal (and its 
aliases) being the orbital period. 

The separation of the 83.38 min signal and the orbital 
period of 81.08 min is conspicuously close to those ob- 
served for the superhu mp period in a la rge number of short 
period dwarf nova ('e.g. |Pattersonl200lj) . In fact, the period 
excess found in HS2331, e = (P sh - P OT b)/Porb = 0.028, is 
similar to that observed in LL And, which has an orbital 
period close to that of HS2331. Based on the observational 
analogies between HS2331 and a number of well-observed 
short-period dwarf novae, it is tempting to ascribe the ob- 
served 83.38 min signal to a permanent superhump. If this 
hypothesis were true, HS2331 would be unique, as it dis- 
plays a permanent quiescent superhump. If 88.51 min were 
the true period its interpretation in terms of a permanent 
superhump would result in e = 0.092, much larger than 
any period excess observed in short period CVs, and we 
might as well have discovered a hitherto unknown phe- 
nomenon. 

3.3. Higher frequencies: White dwarf pulsations &. 
spin ? 

As mentioned above, the periodogram of HS2331 contains 
a number of distinct signals at high frequencies (Fig.EJ), 
most pronounced at 5.61 min (256.6 d^ 1 ) and 1.12 min 
(1284. 7d _1 ). We have inspected in more detail the pe- 
riodogram computed from the 2003 Kryoneri data, which 
provides the longest coverage at a high time resolution, 
as well as from the ULTRACAM observations, which has 
the highest time resolution but covers only ~ 2.5 h. Tabled 
lists the frequencies and amplitudes of all signals that are 
unambiguously detected in the two data sets. The fact 
that most signals are detected over long periods of time 
rules out the CV-typical flickering to explain their origin. 

The extremely complex power spectrum in the range 
around ~ 5 min is indicative of a combination of many fre- 
quencies. Such a structure is reminiscent of the patterns 
detected in the power spectra of ZZ Ceti stars, non-radially 
pulsating white dwarfs with a hydrogen-rich atmosphere. 
They are found in a narrow range of effective tempera- 
tures, 12 500 K > r off > 10 700 K, and show multiperiodic 
variations with amplitudes reaching up t o —0.3 mag an d 
periods ranging from -100 to -1200 s <|Clemensl ll993K 
Resolving the rich frequency spectrum of ZZ Ceti stars re- 
lated to the various pulsation modes, their harmonics, and 
various linear combinations of modes requires alias-free 
uninterrupted mu lti-longitude observing campaigns (e.g. 
iKepler et af]l2003^ . In the case of HS2331, the sampling 
of our 2003 Kryoneri run introduces too much of an alias 
structure (Fig.EJ) to allow a secure identification of all the 
various frequencies. 

Whereas our short ULTRACAM run is not very use- 
ful in identifying the frequencies present in HS2331, these 
data show that the power of the — 5 min variability in- 



Table 3. High-frequencies signals detected in the photom- 
etry of HS 2331+3905. 
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creases towards the blue, which is what would be ex- 
pected for ZZ Ceti pulsations. Untypical for ZZ Ceti pulsa- 
tions are the additional large-amplitude signals detected 
at extremely high frequencies (Fig.EJ). The shortest pe- 
riod observed in any ZZ Ceti star is 70.9 s in G 185-32 
l|Castanheira et aT1l2004|) . Our detection of power at the 
first harmonic of the 1.12 min signal seems to rule out a 
conventional ZZ Cet pulsation as the origin of the high- 
frequency signals. Yet another stable clock available in a 
CV is the spin of the white dwarf, and we tentatively sug- 
gest that the 1.12 min and 0.56 min signals are related to 
the white dwarf spin period and its first harmonic, re- 
spectively. Intriguing is the doublet-like structure of the 
1.12 min and 0.56 min signals (Fig.[f)|). Fast rotation of the 
white dwarf will result in extremely non-linear perturba- 
tions of the pulsation frequencies, and may be responsi- 
ble for this structure. Alternatively, the doublet-structure 
may reflect a sideband frequency, as typically observed 
in intermediate polars. The choice of the lower frequency 
alias in the doublets is somewhat ambiguous: assuming 
that 1277. 7d~ 1 and 2555.5 d^ 1 are the true frequencies 
of the doublet and its first harmonic, the implied beat 
frequency is 7d _1 , which is nearly identical to the spec- 
troscopic period found from the radial velocity variations 
of the Balmer emission line wings (Sect.0J. 

4. Analysis: Spectroscopy 

4.1. Radial velocities 

Radial velocities were m easured using the convolu tion 
technique firs t outlined byfechneider fc Younel l)l980|) and 
developed bv lShafterl (|1983). 
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The emission line profile is convolved with two 
Gaussians of equal FWHM (of the order of the wave- 
length resolution) but with opposite signs. The peaks of 
the Gaussians are separated by a distance chosen in order 
to measure as far out as possible into the wings, while still 
retaining an adequate signal-to-noise level. By doing this 
we ensure that the measured radial velocities trace the 
inner part of the accretion disc, where one may hope to 
find an azimuthally symmetric structure that traces the 
orbital motion of the white dwarf. Prior to the actual ra- 
dial velocity measurement the spectra were transformed 
into the heliocentric rest frame. 

4.1.1. The end of a myth... 

Radial velocity variations were measured for each individ- 
ual observing run (Tabled using the strongest Balmer 
lines as well as He I A5875 and He I 4471. In all cases 
the Gaussians separation was set to values between 2600 
kms -1 and 3000 kms -1 . The radial velocity measure- 
ments were then subjected to the Scargle period analysis 
in the MIDAS/TSA context. The results for each individ- 
ual run did not depend on the emission line used for the 
analysis, and we will discuss in the following the radial 
velocity variations measured in H7 (Calar Alto 2000; we 
opted against H/3 because of the potential contamination 
from the nearby He lines) and Ha for all the other runs. 

Figure|5] shows the periodograms and the correspond- 
ing radial velocity curves folded over the highest peaks 
detected in each of the data sets. The bottom panel shows 
the periodogram and folded radial velocity curve obtained 
when combining all data. The Calar Alto 2000 data show 
a smooth sinusoidal radial velocity variation with an am- 
plitude of ~ 220kms _1 and a period of 3.53 h, which is 
not anywhere close to the orbital period determined from 
the eclipses! This makes HS2331, to our knowledge, the 
first CV in which the dominant radial velocity variation 
is not modulated on the orbital period. With HS2331, 
there are now four CVs in which the photometric modu- 
lation and the spectroscopic modulation are significantly 
different; the o ther three being FS Aur. GWLib and Aqrl 
(see Table 2 of lWoudt et~al]l2004h . The origin of the sec- 
ond period in these systems is much of a mystery as the 
~ 3.5 h periodicity in HS2331, but contrary to HS2331, 
in FSAur, GWLib and Aqrl the spectroscopic variation 
is believed to be orbital. Fig. UHl shows the trailed spectra 
of Ha and H/3 obtained from the observing night with the 
largest data set (i.e. Calar Alto 2003, Sep 11; TablejTJl in 
minutes and phase folded with the orbital period respec- 
tively. The 3.45 h period is clearly seen in the wings of the 
Balmer emission lines. 

The analysis of the other runs reveals that the ra- 
dial velocity variation changes on time scales of several 
days. The structure of the periodograms becomes increas- 
ingly complicated with the length of the run, and fold- 
ing the radial velocity measurements over the strongest 
peak (or any other peak) will not provide a smooth ra- 



dial velocity variation (Fig.^J. The strongest signals in 
the periodograms are found at 4.22 h (INT 2002), 3.45 h 
(Calar Alto 2003), 3.42 h (MDM Hiltner 2003) and 3.45 h 
(all data). All runs, with the exception of the INT data, 
give very close period values which appear consistent with 
each other. Close inspection of the bottom panel of Fig-EI 
shows that, when folded over the ~ 3.45 h signal, (the 
highest peak detected in the periodogram calculated from 
combining all the spectroscopic data together) a shift be- 
tween the individual nights appears (Calar Alto 2000/2003 
and MDM 2003 data). In the case of the INT observing 
run, the poor sampling of the data (Table^) may prevent 
us from detecting the true underlying periodicity at this 
epoch and the relatively smooth appearance of the radial 
velocity curve folded over 4.22 h (the highest peak in the 
periodogram) may be coincidental. Statistically, there is 
no strong preference for any of the aliases found in the 
range ~ 4d _1 — 10 d -1 . 

We conclude that a persistent large-amplitude radial 
velocity variation with a period of ^3.5 h is present in 
HS2331, however, this variation is not coherent and its 
period and/or phase drifts on time scales of days. 

4.1.2. ...sanity refound 

The fact that the dominant radial velocity variation in 
HS2331 does not occur on the orbital period is most sur- 
prising, and alarming, as radial velocity variations mea- 
sured from the emission lines are routinely used not only 
to determine the orbital periods of CVs, but also to es- 
timate the radial velocity of the white dwarf. In order 
to probe for lower amplitude radial velocity variations, we 
have detrended the Calar Alto 2000 and 2003 radial veloc- 
ity measurements with a sine wave adopting the frequency 
of the strongest peak in the corresponding periodogram, 
3.53 h and 3.45 h, respectively. The Scargle periodograms 
of the detrended data (Fig.^1 left panel) contain sig- 
nificant power at frequencies which are consistent with 
the orbital period derived from the photometry (Sect.[3j): 
81.26±0.96min (Calar Alto 2000), 81.08±0.18min (Calar 
Alto 2003), and 81.08 ± 0.04 min (combining both sets). 
Whereas the signal of the orbital period completely dom- 
inates the periodogram of the 2000 data, significant resid- 
ual power is left near 3.45 h in the case of the 2003 data. 
The main difference between the two data sets is the ex- 
tent of the observations. While the 2000 data were ob- 
tained during two consecutive nights, the 2003 data was 
obtained during three nights spanning a total of six nights. 
As the dominant radial velocity variation drifts on times 
scales of a few days, dctrcnding with a sine with a fixed 
period will become increasingly imperfect with increas- 
ing length of the observing run. Hence, the residual low- 
frequency power in the case of the 2003 data reflects the 
imperfection of the applied detrending. 

Folding the detrended radial velocities over the orbital 
period (i.e. 81.08 min) results in a sinusoidal radial velocity 
curve (Fig.^1 right panel), which is smoother for the 2000 
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data than for the 2003 data, probably again a result of 
the cleaner detrending possible for these observations. The 
consistent detection of the 81.08 min period both, in the 
photometry and in the spectral radial velocity variations 
lends further support for our conclusion that this is indeed 
the orbital period of HS 2331. 

A sine fit to the folded detrended radial velocities of 
the Calar Alto 2000 data (top right panel of Fig.[j*l"| re- 
sults in an amplitude of ~ 32kms _1 , a systemic velocity 
of ~ lkms , and a phase offset with respect to mid 
eclipse of ~ 13°. The amplitude, interpreted as radial ve- 
locity of the white dwarf Ki as it is commonly done, is 
entirely c onsistent with what is observ ed in short-period 
CVs ('e.g. iThorstensen fc Fentonll2003h . The phase offset 
with respect to the expected motion of the white dwarf is 
most likely related to the fact that the eclipse is caused by 
an obscuration of the hot spot in the accretion disc, and 
not of the white dwarf itself. In addition, even in "well- 
behaved" CVs the broad wings of the emission lines do not 
necessarily exactly trace orbit of the white dwarf, most 
likely because of slight asymmetries in the disc. 

4.2. A model for the spectral energy distribution 

We have combined our optical data with the 2MASS J, 
H, and K magnitudes of HS2331 (FigJEJ, and mod- 
eled the optical/IR spectral energy distribution (SED) 
with a three-component model comprising: (1) a synthetic 
white dwarf spectrum to account for the observed broad 
Balmer absorption lines. A grid of model spectra cover- 
ing T cff = 8 000 - 15 000 K was calculated with TLU STY 
and SYNSPEC !lHubenvlll988llHubenv fc Lanzlll995h . We 
fixed the surface gravity to logg = 8.0, corresponding to 
M w d ~ 0.6 Mq, and the photospheric abundances to 0.1 
of their solar values. (2) The emission of an isothermal 
and isobaric slab with a finite depth to account for the 
observed Balmer emission lines and any as sociated contin- 
uum. This model is described in d etail in IGansicke et all 
l|l997(l and IGansicke et all l)l999j) . In brief, the free pa- 
rameters for the generation of this component were the 
temperature of slab and the column density of the slab 
along the line of sight. The temperature accounts for the 
ionisation/excitation of the lines and the optically thick 
blackbody envelope, and the column density determines 
the ratio in optical depth between continuum and lines, 
i.e. the equivalent width of the lines. (3) A late-type spec- 
tral template to account for the emission of the secondary 
star in HS2331. The libr ary of templates c o vered spectral 
types M0. 5 to M9 from iBeuermann et all lll998t) and LP 
to L8 from lKirkpatrick et a l. (1999) and Kirkpatrick et alJ 
( 2000j. The L-dwarf type stars were normaliz ed to surface 
fluxes using the radii and distances given bv iDahn et alJ 
l|2002h . the M-dwarf templates provided by Beuermann 
were already normalized in this way. 

We started by approximately modeling the observed 
broad Balmer absorption lines with a white dwarf model, 
which suggested T c g ~ 10 000 K and a white dwarf flux 



contribution of - 80% at 5000 A (Fig-EJ. The next step 
was to extrapolate the white dwarf model established in 
this way into the FUV wavelength range covered by STIS 
(FigJlflp. Interpreting the sudden flux up-turn at wave- 
lengths A <; 1650 A as photospheric emission from the 
white dwarf we refine the white dwarf temperature to 
T cff ~ 10 500 ± 250 K. The red end of the STIS range 
is very sensitive to the white dwarf temperature because 
of the strong temperature dependence of the quasimolcc- 
ular H2 absorption at 1600 A. The flux scaling factor of 
our white dwarf model implies a distance of ~ 90 pc. The 
systematic error in this distance estimate depends primar- 
ily on the assumed white dwarf mass, assuming a mass in 
the range 0.4 - 0.9 M Q results in d = 90 ± 15 pc. A low 
distance to HS2331 is supported by the large proper mo- 
tion of the star, fj, ~ 0.14"/yr _1 , which has been derived 
from the comparison of the DSS1 and DSS2 images of the 
system. An isothermal/isobaric slab with a temperature of 
T = 6500 K and a surface density of £ = 1.8 X 10 -2 gcm~ 2 
accounts well for the observed Balmer emission and the 
Balmer jump, and contributes ~ 20% of the flux at 
5000 A. Assuming d — 90 pc the area of this simple "disc" 
component is comparable to the size of the white dwarf 
Roche lobe. The red end of our optical spectrum shows no 
evidence for the TiO absorption bands that are observed 
in CVs where the late-type donor star significantly con- 
tributes to the red/near-IR spectrum. In the last step of 
our SED modeling we add the spectral contribution of the 
secondary, assuming R2 — 8 x 10 9 cm (corresponding to 
-Porb = 81.08 min and q = M w d/M scc = 9; see Sect. l5.1JI . 
The flux limit imposed by the 2MASS J, H, and K mea- 
surements constrain the spectral type of the secondary to 
be ~L2 or later (Fig.[HJ. 

Whereas our three-component model of the FUV-to- 
IR SED of HS2331 is overall rather satisfying, two re- 
gions of excess flux are observed: the FUV continuum be- 
low ~ 1600 A significantly exceeds the predicted white 
dwarf flux (the simple "disc" model and the secondary 
contribute no flux in the FUV), and the optical contin- 
uum in the range ~ 7000 - 9000 A exceeds the flux of the 
sum of all three components. We believe that the most 
likely explanation for this excess flux is that our model 
for the disc emission is too naive - we just assume a single 
temperature/density pure hydrogen slab. The FUV emis- 
sion lines clearly indicate that the disc contains hotter 
regions, which will contribute to some extent to the FUV 
continuum. The excess flux observed at the red end of our 
optical spectrum may be related to the Paschen contin- 
uum of a somewhat colder part of the disc. It is likely 
that a multi-temperature/density model of the quiescent 
disc could account for the observed excess flux, however, 
we feel that the number of additional free parameters does 
not warrant this exercise. 

5. Discussion 

Before discussing the nature of the system, we briefly sum- 
marise the results of our analysis. 
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Fig. 13. Gray line: the STIS FUV spectrum of HS2331 
from Fig.03 Black line: the white dwarf model from Fie . 1 121 
extended into the STIS range. The spectrum of a ~ 
10 500 K white dwarf reproduces well the sudden increase 
of flux observed longwards of 1600 A. 



1. HS2331 has been identified as a CV because of its emis- 
sion lines detected in the HQS. 

2. Coherent low-amplitude eclipses detected over four 
years unequivocally determine the orbital period to be 
P orb = 81.08 min. The mean light curve of HS2331 
bears extreme similarity to that of WZ Sge. 

3. The optical and FUV spectrum of HS2331 is domi- 
nated by the photospheric emission of a white dwarf 
with T cff ~ 10 500 K. 

4. The power spectrum of our optical photometry con- 
tains a number of significant signals in the range 
5 — 6 min, which we interpret as non-radial ZZ Ceti 
pulsations. 

5. Additional signals are found at periods of 1.12 min and 
0.56 min, which we tentatively identify as the spin pe- 
riod of the white dwarf and its first harmonic. 

6. The distance of the system is d ~ 90 ± 15 pc, derived 
from modelling the emission of the white dwarf. This 
low distance is consistent with the high proper motion 
of the system, /i ~ 0.14"/yr~l. 

7. Combining our optical spectroscopy of HS2331 with 
the 2MASS J, H, and K magnitudes of the sys- 
tem, and adopting d — 90 pc and a mass ratio of 
q = M w( j/M sec = 9 (see below, Sect. l5.lfl . we find that 
the donor star must be of extreme late spectral type, 
L2 or later. 

8. Our extensive optical photometry consistently con- 
tains a strong signal at a period somewhat longer than 
P or b ■ Selecting the alias of that pattern which is closest 
to P or b gives 83.38 min, which is 2.8 % longer than the 
orbital period. 

9. Our spectroscopic data set reveals consistently the 
presence of large-amplitude radial velocity variation 
with a period of ~ 3.5 h. However, this radial veloc- 



ity modulation is not strictly coherent, but appears to 
drift in either period or phase on a time scale of days. 
10. After detrending the radial velocity measurements 
with a 3.5 h sinusoid, a radial velocity variation with 
a period of 81.08 ± 0.04, i.e. P or b, is detected with an 
amplitude of ~ 32kms~ 1 . 

5. 1 . The evolutionary state of HS 2331 +3905 

A number of arguments suggest that HS2331 is a CV in a 
late stage of its evolution. Its orbital period is close to the 
observed period minimum of CVs. A very low mass trans- 
fer rate is indicated by the low temperature of the white 
dwar f, which is a measure of th e secular mean accretion 
rate l|Townslev fc BildsteiJ l20031 . This is consistent with 
the absence of observed outbursts. All our detailed pho- 
tometric and spectroscopic campaigns carried out since 
September 2000 have found HS2331 near V ~ 16.5. The 
CCD-equipped 0.37 m Rigel Telescope of the University 
of Iowa has been used to monitor HS2331 on 16 nights in 
October to December 2003 and on 15 nights in May/ June 
2004, finding the system at an instrumental magnitude of 
16.4 ± 0.1. Additional visual monitoring of HS2331 has 
been carried out during ~ 250 nights since March 2003 
using a 8" Schmidt-Cassegrain telescope, with no gaps in 
the coverage being longer than one week. On no occasion 
has the star been detected, with a limiting magnitude of 
~ 13.5 ±0.5. In many aspects, HS2331 resembles the well- 
studied short-period dwarf nova WZ Sge, which shows ex- 
tremely bright superoutbursts every 2-3 decades. If this 
analogy were correct, we would predict an outburst mag- 
nitude for HS2331 of V ~ 9 - 10. 

Yet another indication of HS2331 being an old evolved 
CV is the very late spectral type of the secondary star. 
Assuming a white dwarf mass of 0.6 Mq and a secondary 
mass of 0.07 Mq (as it may be realistic for M scc near the 
minimum period), i.e. q ~ 9, and a distance of 90 pc, our 
model for the SED of HS2331 suggests a spectral type later 
than L2 fScct. l4~2*)l . which makes this system one of the 
best candidates for a brown dwarf donor. Detailed infrared 
spectroscopy will be required to test this hypothesis. 

Using the same stellar mass assumption and an incli- 
nation of ~ 75°, based on the detection of grazing eclipses 
in the optical light curve, the radial velocity variation 
of the white dwarf is expected to be K\ ~ 36kms _1 , 
which is consistent with the amplitude of the 81.08 min 
radial velocity variation detected in the spectroscopy, af- 
ter detrending the data for the dominant 3.5 h modulation 
(Sect.EESJ. 

5.2. A permanent superhumper? 

From our four years of photometric data, HS 2331 ap- 
pears to be a permanent superhumper with a likely su- 
pcrhump period of P srl = 83.38 min, i.e. a period excess 
(P s h — P or b)/Porb = 2.8%. This finding is exceptional, as 
superhumps are normally detected during superoutbursts 
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of SU UMa type dwarf novae or in some nova-like vari- 
ables during their high states, i.e. always during phases of 
high mass transfer in the accretion disc. Superhumps are 
thought to be related to the precession of tidally unsta- 
ble eccentric acc retion discs in close binarie s with extreme 
mass ratios (e.g. IOsakilll985t lMurravlll998|) . 

If our interpretation of the photometric signal is 
correct, and indeed P s h = 83.38 min, the i mplied pe- 
riod excess of 2.8% is not particularly small I Pattersor] 
l200lh . somewhat arguing against an extreme mass ra- 
tio in HS2331. However, considering the unique nature 
of HS2331, it is not clear if Patterson's relation can be 
directly adopted for this system. 

The strangest and most disconcerting observational 
fact about HS2331 is without doubt the detection of a 
dominant radial velocity variability with a period of 3.5 h, 
which does not correspond to the orbital motion of the 
system. This periodicity resulted from measuring changes 
in the wings of the Balmcr emission lines and therefore 
reflects variability in the inner disc region of the sys- 
tems, however we have at present no explanation for the 
physical origin/significance of its 3.5 h period. Again, in a 
very speculative approach, we suggest that the period (or 
phase) drift of this radial velocity variation may be related 
to the presence of permanent superhumps. The precession 
period of the accretion disc implied by the orbital period 
and the superhump period is ~ 2 d, which is comparable 
to the time scale on which the 3.5 h radial velocity varia- 
tion varies. An intensive spectroscopic campaign covering 
several of the putative precession cycles could test this 
hypothesis. 

Another important conclusion from the radial veloc- 
ity analysis is that care should be taken when implicitely 
assuming that the dominant radial velocity variation de- 
tected in cataclysmic variables (and related objects) re- 
flects the orbital period. Whereas this is contrary to the 
customary approach for the past decades, the results ob- 
tained for HS 2331 are a clear warning that this custom 
could occasionally lead to wrong interpretations. 

5.3. The brightest CV white dwarf pulsator 

Until recently, only a single CV containing a white dwarf 
pulsa t or ha s been discovered: GWLib ijvan Zvl et alJ 
l200(l 120041) . GWLib is only partially understood, due 
to the lack of a proper mode identification of the ob- 
served pulsation frequencies, and more severely, due 
to the circumstance that the temperature of its white 
dwarf, T e g = 14 700 K, is actually significa ntly above the 
blue end of the ZZ Ceti instability strip ijSzkodv et alJ 
2002J). Three additional CVs containing pulsating ac- 
creting white dwarfs have been discovered in the SDSS: 
SDSS J013132.39-090122.3, SP SS .1161033.64-010223.3 
and SDSS J220553. 98+ 115553. 7 tWarner fc Woudtll2003t 
IWoudt fc Warnerll2004h . 

We have shown that HS2331 exhibits the photometric 
behaviour typical of ZZ Ceti pulsators, showing a multi- 



periodic and complex power spectrum in the range ~ 60 s 
to - 300 s. With V ~ 16.5, HS2331 is by far the brightest 
CV white dwarf pulsator, making it the prime target for 
detailed follow-up studies. The white dwarf temperature 
implied by the SED of HS2331 (T off ~ 10 500 K, assuming 
M w( j — 0.6 Mq, Sect. l4.2[l is just slightly below the red end 
of the ZZ Ceti instab ility strip for single white dwarfs (e.g. 
iBereeron et al.ll2004jh The ultra-short periods detected in 
HS2331, with periods of 1.12 min and 0.56 min which are 
too small for conventional (/-mode pulsations, remains a 
puzzle. One possible solution is to assume that 1.12 min 
is the white dwarf spin period, and 0.56 min its first har- 
monic. Polarimetry and X-ray detection of the 1.12 min 
period would be the ultimately proof of this hypothesis. 

A more detailed understanding of the white dwarf 
pulsations in CVs is just at its beginning, as the first 
models for the drivin g mechanism are being developed 
l|Townslev et al.ll2004j) . and a larger sample of these ob- 
jects is discovered. Future potential asteroseismological 
applications in CV white dwarfs include measuring the 
white dwarf masses, as well as studying the accreted en- 
velope. 

6. Conclusion 

We have discovered a new short orbital period CV, 
HS2331, which displays an overwhelming wealth of ob- 
servable phenomena, and we encourage detailed follow-up 
studies over all wavelength ranges and time scales to ex- 
ploit the opportunities that this system offers for an im- 
proved understanding of CV evolution and accretion onto 
white dwarfs. 
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Fig. 6. Scargle periodograms computed from the CCD photometry of HS2331. Top panel: The entire frequency 
range sampled by the Kryoneri 2003 data, which represents the largest and best sampled data set to investigate 
periodicities on time scales of ~ 2d to ~ 30s. Middle panels: Close-ups of the power spectra calculated from all 
available photometric data (Tabled]). The left panel shows a coherent signal at 35.52 d _1 . Shown in gray is the window 
function of the entire data set, shifted to the frequency of the strongest peak. We interpret this signal as twice the 
orbital frequency, as the orbital light curve is dominated by a double- hump structure (Fig.[7|). The right panel shows 
several signals in the low-frequency range. Practically no power is detected at the orbital frequency (i.e. 17.76 d _1 
or 81.08 min), but a strong signal is found at 17.27d _1 (83.38 min) and the one-day alias at 16.27d _1 (88.51 min). 
Bottom panels: Close-up views of the high-frequency signals in the Scargle periodogram computed from the Kryoneri 
2003 data only. As in the middle panels, the gray line shows the alias structure resulting from the sampling of our 
observations. Clear signals are detected at 5.61min (256.6d _1 ), 1.12min (1284. 7d _1 ), and 0.56 min (2555. 5d _1 ). 
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Fig. 7. HS2331 entire photometric data set folded over Fig. 8. HS2331 entire photometric data set folded over 
81.08 min (upper panel) and re-sampled into bins of 100 83.38 min (upper panel) and re-sampled into bins of 100 
points (bottom panel). The light curves are repeated over points (bottom panel). The light curves are repeated over 



two cycles for clarity. 



two cycles for clarity. 
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Fig. 9. Left Panel: Scargle periodograms of the individual data sets (see TableQJ and of all the data combined together. 
Right Panel: Radial velocity curves of the individual data sets and of all the data sets combined together, folded over 
the highest peak in their respective Scargle periodograms (shown on the left). The zero phase is taken to be the first 
point in each of the data sets respectively and the curves are repeated twice for clarity. 
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Fig. 10. From left to right: Trail spectra of Ha and H/3 computed from the largest data set within an observing night 
(i.e. Calar Alto 2003, Sep 24; see Tabled, and the same Ha and H/3 trail spectra but this time folded over the orbital 
period (i.e. P or b = 81.08 min). The white stripes in the trail spectra of the left panel are caused by lack of observations 
during the corresponding times. 
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Fig. 11. Left Panel: From top to bottom, Scargle periodograms of the detrended Calar Alto 2000 and Calar Alto 
2003 spectroscopic data sets and a combination of both sets (i.e. Calar Alto 2000 & Calar Alto 2003). The radial 
velocity measurements were detrended by subtracting a sine wave with their respective highest peak value (see the 
left panel of Fig.|3) as the frequency. The combined Calar Alto 2000 & 2003 Scargle periodogram in the bottom panel 
was computed after detrending the data sets individually. Right panel: Detrended radial velocity curves folded over 
the orbital period obtained from the photometry in Sec. |3] (i.e. 17.76 d" 1 or 81.08 min). The zero phase in these curves 
corresponds to the photometric zero phase as defined in Sec. Inland the curves are repeated twice for clarity. Also shown 
on the top right panel is the best fit to the detrended radial velocity curve (see text for more details). 




Fig. 12. Three component model to the optical-IR spectral energy distribution of HS2331. Plotted in gray are the 
observed average spectrum from Fig.0]along with the 2MASS J, H, and K fluxes (filled circles). Plotted in dotted lines 
are the white dwarf model (open squares denote the IR fluxes), the isothermal slab model representing the optically 
thin disc emission (stars denote the IR fluxes), and the spectral contribution of a L2 donor star (open circles denote 
the IR fluxes). The sum of the three components is shown as black solid line and as black filled circles. 



